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Preface 



Differential equations of the fractional order appear more and more frequently 
in different research areas and engineering applications. An effective and easy- 
to-use method for solving such equations is needed. 

However, known methods have certain disadvantages. Methods, described 
in details in ^, [| for fractional differential equations of the rational order, do 
not work in the case of an arbitrary real order. On the other hand, there is an 
iteration method described in ||, which allows solution of fractional differential 
equations of an arbitrary real order, but it works effectively only for relatively 
simple equations, as well as the series method [l], pTj . Other authors (e.g. 
ID, ||) used in their investigations the one-parameter Mittag-Leffler function 
E a (z) = J2k?=o TJak+Tj- Still other authors [||, fj]] prefer the Fox //-function ||, 
which seems to be too general to be frequently used in applications. 

Instead of this variety of different methods, we introduce a method which is 
free of the mentioned disadvantages and suitable for a wide class of initial value 
problems for fractional differential equations. The method uses the Laplace 
transform technique and is based on the formula of the Laplace transform of 
the Mittag-Leffler function in two parameters E a g(z). We hope that the de- 
scribed method could be useful for obtaining solutions of different applied prob- 
lems, appearing in physics, chemistry, electrochemistry, engineering, financing 
and banking, etc. To outline the area of the method's applicability, we have 



included in the bibliography also the works by different authors |30|| — [ 54 1 , in 
which fractional linear differential equations appears or which could serve as a 
basis for obtaining such equations. 

This work deals with solution of the fractional linear differential equations 
with constant coefficients and consists of four short chapters. 

In Chapter 1 we present some auxiliary tools which are necessary for us- 
ing the method. The reader can find the definition there and some important 
properties of the Mittag-Leffler function in two parameters and the Wright func- 
tion. The basic result, presented in Chapter 1, is the Laplace transform of the 
Mittag-Leffler function and its derivatives. Besides that, we introduce two tools 
necessary for testing candidate solutions by direct substitution in corresponding 
equations: fractional derivatives of the Mittag-Leffler function and a rule for 
the fractional differentiation of integrals depending on a parameter. 

In Chapter 2 we give solutions to some initial- value problems for "standard" 
fractional differential equations. Some of them were solved by other authors 
earlier by other methods, and the comparison in such cases just underline the 
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simplicity and the power of our approach. 

In Chapter 3 we extend the proposed method for the case of so-called "se- 
quential" fractional differential equations (we adopted the convenient terminol- 
ogy of Miller and Ross). For this purpose, we obtained the Laplace transform 
for the "sequential" fractional derivative. The "sequential" analogues of the 
problems, solved in Chapter 2, are considered. Naturally, we arrive at solutions 
which are different from those obtained in Chapter 2. 

However, there is something common in solutions of the corresponding 
"standard" and "sequential" fractional differential equations: they both have 
the same fractional Green's function. In Chapter 4 we give our definition of the 
fractional Green's function and some of its properties, necessary for construct- 
ing solutions of initial- value problems for fractional linear differential equations 
with constant coefficients. 

We give the solution of the initial-value problem for the ordinary fractional 
linear differential equation with constant coefficients using only its Green's func- 
tion. Due to this result, the solution of such initial value problems reduces to 
finding the fractional Green's functions. We obtained the explicit expressions 
for the fractional Green's function for some special cases (one- , two- , three- 
and four-term equations). 

The explicit expression for an arbitrary fractional linear ordinary differential 
equation with constant coefficients ends this work. 

In what follows below, a Df means the Riemann-Liouville fractional deriva- 
tive of order a (e.g., Q, 

The author thanks Ms. Serena Yeo for checking the language of this work. 



4 



Chapter 1 

Introduction 



1.1 Mittag-Leffler function in two parameters 

The function, which plays a very important role in this work, was in fact in- 
troduced by Agarwal [23]. A number of relationships for this function were 
obtained by Agarwal and Humbert [^4|] using the Laplace transform technique. 
This function could have been called the Agarwal function. However, Agar- 
wal and Humbert generously left the same notation as for the one-parameter 
Mittag-Leffler function, and that was the reason that now the two-parameter 
function is called the Mittag-Leffler function. We will use the name and the 
notation used in the fundamental handbook on special functions [1C]. In spite 
of using the same notation as Agarwal, the definition given there differs from 
Agarwal's definition by a non-constant factor. 

Definition. A two-parameter function of the Mittag-Leffler type is defined by 
the series expansion [10] 

°° k 

^ w -£rtS+fl- (a>0 ' " >0 > (u > 



Relation to other functions. There are some relationships given in [10, 
§18.1]: 

E 1:1 (z)=e z , E lt2 {z) = e -^, (1.2) 

E 2>1 (z) = cosh(^), E 2>2 (z) = Smh y i) , (1.3) 

E 1/2 ^z) = -^=e-*erfc(-v^). (1.4) 



For (5 = 1 we obtain the Mittag-Leffler function in one parameter: 

oo k 

Wz)= SrfrST!) sE " (2) ' (L5) 
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The function £t(u,a), introduced in Hj for solving differential equations of the 
rational order, is the particular case of the Mittag-Leffler function 

00 (nt) k 
^0 F ^ + k + X ) 



Rabotnov's (TJ| function 3 a (/3,t) is the particular case of the Mittag-Leffler 
function ( |1.1|) too: 

00 ak.Lk(a+l) 

B a (M = t Q E r(( / +1)(1 + a)) = (1-7) 

It follows from the relationships Ql.6| ) and ( |1.7[ ) that the properties of the 
Miller-Ross function and Rabotnov's function follows from the properties of the 
Mittag-Leffler function in two parameters (|1.1|). 



Plotnikov cf. [26|] and Tseytlin |g6[ used in their investigations two 
functions Sc a (z) and Cs a (z), which they call the fractional sine and cosine. 
Those functions are also just the particular cases of the Mittag-Leffler function 
in two parameters: 

00 C_1 )nj2-a)n+l 
00 (_1 \n„(2-a)n 

Of course, properties of the fractional sine and cosine follow from the prop- 
erties of the Mittag-lefner function (|l]J) . 

1.2 The Laplace transform of the 

Mittag-Leffler function in two parameters 



As follows from relationship (1.2), the Mittag-Leffler function E a n(z) is a gen- 
eralization of the exponential function e z and the exponential function is a 
particular case of the Mittag-Leffler function. 

We will outline here the way to obtain the Laplace transform of the Mittag- 
Leffler function with the help of the analogy between this function and the 
function e z . For this purpose, let us obtain the Laplace transform of the function 
t k e at in an untraditional way. 

First, let us prove that 



00 1 

t+zt-- 1 



e- l e ±zt dt = , \z\ < 1. (1.10) 

1 T z 



Indeed, using the series expansion for e 2 , we obtain 



1 ^ i±z 
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OO 



e-tgtrfi = = V / e-H k dt = V(±z) fc = -— - — . (1.11) 
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Second, we differentiate both sides of equation ( 1.10 ) with respect to z. The 
result is 

l e-H^dt = — - FTT , (|z|<l), (1.12) 



and after obvious substitutions we obtain the well-known pair of the Laplace 
transform of the function t k e *: 



poo 

/ e- pt t k e ±at dt 
Jo 



kl 



,fc+l 



(Re(p) > \a\). 



(1.13) 



Let us now consider the Mittag-Lefner function (|l.l| ). Substitution of ( |l.l| ) 
in the integral below leads to 



roo 1 

/ e-HP- x E a , p {zt a )dt = - , 

Jo 1 — z 



M <i), 



(1.14) 



and we obtain from (|1.14| ) a pair of the Laplace transforms of the function 



(fc). 



A/ 



ft! 



e - P t t ak+P-l E (k) r± at aj dt ^ 



(i2e(p) > |a| 1/Q ). (1.15) 



The particular case of (|1.15| ) for a = (5 = ^ 

fc! 



f°° e~P t t ! ^ 1 E^ 1 (±aVi)dt 

JO 2>2 



is useful for solving semidifferential equations considered in M, 



(fle(p) > a 2 ). (1.16) 



1.3 The Wright function 

This function, related to the Mittag-Lefner function in two parameters E a ^(z), 
was introduced by Wright p9| , cf. Jill 24]]. A number of useful relationships was 
obtained by Humbert and Agarwal [24] with the help of the Laplace transform. 
For convenience we adopt here Mainardi's notation for the Wright function 

W(z;a,/3). 

Definition. The Wright function is defined as [10, formula 18.1(27)] 



k=0 



(ak + (3) 



(1.17) 



Integral representation. This function can be represented by the following 



integral |Tg formula 18.1(29)]: 

W(z\a,P) = — f u- p e u+zu ~ a du (1.18) 
2iri J Ha 

where (Ha) denotes Hankel's contour. 
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To prove ( |1.18| ), let us write the integrated function in the form of a power 
series in z and perform term-by-term integration using the well-known formula 
§ 1-6(2)]: 

i 

e u u z du. 



T{z) 



Ha 



Relation to other functions. It follows from the definition ( 1.17 ) that 

W(*;0,1) = e z 



- W(t-; 1,^ + 1) 



Ju(z) 



Taking (3 = 1 — a, we obtain Mainardi's function M(z] a): 

_ a , 1 - a ) = M(z; a) = £ ^rTTTT^ 

^fc! r(-a(fc + l) + 1) 



(1.19) 
(1.20) 

(1.21) 



The following particular case of the Wright function was given by Mainardi 



W(-z; 



1 1. 



M(z; 



IT 



exp 



(1.22) 



We see that the Wright function is a generalization of the exponential func- 
tion and the Bessel functions. For a > and (3 > it is an entire function in z 
©• 

Recently Mainardi (T^j pointed out that W(z;a, (3) is an entire function in 
z also for —1 < a < 0. 

Let us prove this statement. Using the well-known relationship 

7T 



T(y)T(l-y) 



sin(-7T7/) ' 

we can write the Wright function in the form 



W(z;a,(3) = -J2 
vr f — „ 



1 ~ z fcr(l -ak- P) sinvr(aA; + /?) 



fc=0 



A;! 



Let us introduce an auxilary majorizing series 



1 OO 

s = -y 



r(l -ak-0) 



The convergence radius of series ( |1.25| ) for — 1 < a < is infinite: 



R = lim 

fc^oo 



T(l-ak-(3) (k + l)\ 



k\ 



k + 1 
urn | — j — —— 

fc^oo a \ a k c 



oo. 



(1.23) 



(1.24) 



(1.25) 



(1.26) 



r(l — ak — a — f3) 

(We use here relationship || formula 1.18(4)].) 

It follows the comparison of the series ( p7D and (|L25| ) that for a > -1 and 
arbitrary /? the convergence radius of the series representation of the Wright 
function W(z;a, f3) is infinite, and the Wright function is an entire function. 

The Wright function plays an important role in the solution of linear partial 
fractional differential equations, e.g. fractional diffusion-wave equation. 
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1.4 Tools for testing candidate solutions 



The following tools are necessary for checking by substitution if candidate so- 
lutions, obtained by any more or less formal method, satisfy corresponding 
equations and initial conditions. 

Fractional derivatives of the Mittag-Lemer function. 

By fractional-order differentiation qD^ (7 is an arbitrary real number) of series 



representation (1.1) we obtain 



oD 7Q*k+f>-l E m ( At «)) = t ^+/3-7-l£( fc )_ 7 ( A r) (1.27) 



The particular case of relationship (1.27) for k = and integer 7 is given in 



|1C[, equation 18.1(25). 



A rule for fractional differentiation of an integral depending on 
a parameter, when the upper limit also depends on the parameter: 







Df I K{t,r)dT = f T D?K(t,T)dT + lim T D^ l K{t,r), (0 < a < 1). 
Jo Jo 

(1.28) 

The proof is straightforward: 



I d f x [ x K(t,0dt 



T(l — a) dx Jo J£ {x — t) a 

d r x ~, ^ r x d 



/ K(x,0d£= / —K{x,£)d£+ lim K(x,£) 
Jo Jo dx z^x-o 

[\d«K(x, e)de + lim zD«~ l K(x, 0, 
Jo i^x-o 



dx 
rx 



where 



K(x,0 



x K(t,i)dt 



T(l - a) Ji (x - t) a ' 
The following important particular case must be mentioned. If we have 



K(t — t)/(£) instead of K (t, r), then relationship ( |1.28| ) takes the form: 

B? fK{t - r)f(r)dr = f^K{j)f{t - r)dr + lim f(t - r) Q D^ l K(r)) 
Jo Jo 

(1.29) 
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Chapter 2 

Standard fractional 
differential equations 



The following examples illustrate the use of ( 1.15| ) for solving fractional-order 
differential equations with constant coefficients. In this chapter we use the 
classic formula for the Laplace transform of the fractional derivative, as given, 
e.g., in Q, p.134 or §, p.123: 



/ e~* D?f(t) dt = p a F(p) - £ p k D 



t=o 



(n 



1 < a < n). 
(2.1) 



2.1 Ordinary linear fractional 
differential equations 

In this section we give some examples of solution of ordinary linear differential 
equations of the fractional order. 



Example 1. A slight generalization of an equation solved in (|l|], p. 157): 



Dl /2 f(t) + af(t) = 0, (t>0); 

Applying the Laplace transform we obtain 

C 



D; 1/2 f(t) 



t=o 



C 



(2.2) 



F{p) 



c= D; 1/2 f(t) 



t=o 



p 1 / 2 + a 

and the inverse transform with a help of ( |1.16| ) gives the solution of (|2.2|): 

f(t) = Cr l / 2 E L i(-oVt). (2.3) 



Using series expansion ( |l.l|) of E a> p(t), it is easy to check that for a = 1 solution 
( f2li| ) is identical to the solution f(t) = C(-^ — e*erfc(\/i)) , obtained in Q in 



(2.4) 



a more complex way. 

Example 2. Let us consider the following equation: 

D?f(t) + Dff(t) = h(t), 
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which "encounters very great difficulties except when the difference q — Q is 
integer or half-integer" (Q, p. 156). 

Suppose that < q < Q < 1. Laplace transform of fl2.4|) leads to 



(pQ +p i)F(p) = C + H{p), where C = (oA^V(i) + oA 



Q-l 



i=0 



and 



F(p) 



C + ff(p) _ C + ff(p) 



(C + tf(p)) 



(2.5) 



(2.6) 



pQ+p<i pi(pQ-i + l) v ^" pQ-i + 1 

After inversion with a help of (|1.15) for a = Q — q and /? = Q, we obtain 



the solution: 



/(t) = C G(t) + f G(t - r)h(r)dr, 
Jo 

c = (ourvw + oa q_1 /w) | t=0 , g(o = 



(2.7) 



The case < q < Q < n (for example, equation obtained in p2|j) can be 
solved similarly. 

Example 3. Let us consider the following initial value problem for a non- 
homogeneous fractional differential equation under non-zero initial conditions: 



B?y(t) - Xy(t) = h(t), (t > 0); (,£>?"*!/(<) 



t=o 



b k , (fe = 1,2, 



,n 



(2.8) 

where n — 1 < a < ra. Problem fl2.8| ) was analytically solved in || by the 
iteration method. With the help of the Laplace transform and formula ( 1.15| ) 
we obtain the same solution directly and more quickly. 

Indeed, taking into account initial conditions ( |2.8D , the Laplace transform 
of equation (ET1) yields 



P a Y( P ) - XY(p) = H(p) + hp 

k=l 



k-1 



Y{p) 



H{p) i p k - x 



(2.9) 



and the inverse transform using (1.15) gives the solution: 

n ft 

yit) = Y,bkt a - k E ata _ k+ i(\t a )+ / (t-r) a - l E a)a (\(t-T) a )h{T)dT. (2.10) 

Jo 



k=l 



11 



2.2 Partial linear fractional 
differential equations 

The proposed approach can be successfully used for solving partial linear dif- 
ferential equations of the fractional order. 

Example 4. Nigmatullin's fractional diffusion equation. 

Let us consider the following initial boundary value problem for the frac- 
tional diffusion equation in one space dimension: 



D?u(x,t)=\ 



,d 2 u(x,t) 
dx 2 



(t>0, 



lim u(x,t)=0; D^~ l u{x,t) 

x— >±oo 



oo < x < oo); 

t=0 



(2.11) 
(2.12) 



We assume here < a < 1. Equation of the type (|2,ll| ) was deduced by 
Nigmatullin [ 13 1 and by Westerlund [14] and studied by Mainardi |jl5|1). We 
will give a simple solution of problem (2.11) demonstrating once again the 
advantage of using the Mittag-Lefner function in two parameters ( |1 . 1| ) . 

Taking into account boundary conditions ( |2.12 ), the Fourier transform with 
respect to variable x gives: 



D?u(f3, t) + X 2 p 2 u(P, t)=0 



ia- 1 



D?- L u(x,t) 



t=0 



HP), 



(2.13) 
(2.14) 



where (3 is the Fourier transform parameter. Applying Laplace transform to 
( p. 13 ) and using initial condition ( |2.14| ) we obtain 



U(J3,p) 



p a + A 2 /? 2 ' 

Inverse Laplace transform of ( 2.15| ) using ( 1.15| ) gives 
u{(3,t) = ip{(3)t a - 1 E ol J-\ 2 2 t a ), 



(2.15) 



(2.16) 



and inverse Fourier transform produces the solution of the problem ( 2.11[ )- 



/oo 
G{x-C,tMC)dC, 
-oo 



G(x,t) 



1 

7T Jo 



t a ~ E aa (^ — A p t") cos f3xd/3 



(2.17) 
(2.18) 



Let us evaluate integral ( 2.18| ). The Laplace transform of ( p.18 ) and formula 
1.2(11) from @ produce 



g(x,p) 



1 



cos(Px d(5) 



IT J A 2 /? 2 + p 
and the Laplace inversion gives 

1 



J_ -a/2 |x|A-V /2 

2A P 



G(x,t) 



AXtti 



e pt p 2 exp(-|x|A 1 p a ^ 2 )dp. 



(2.19) 



(2.20) 



Br 
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Performing substitutions a = pt and z = \x\X~ 1 t~ p (p = a/2) and trans- 
formation of the Bromwich contour (Br) to the Hankel contour (Ha), as was 
done in a similar case by Mainardi |l5|], we obtain 

c <* t >-£ss/^£-sr , " r <-*-»'>- * = S (2 ' 21) 

_ffa 



where W(z,X,p) is the Wright function ( |1,17 ). We would like to note that, 



in fact, we have just evaluated the Fourier cosine-transform of the function 
Ul ((3) =t a - 1 E a , a (-\ 2 p 2 t a ). 

It is easy to check that for a = 1 (traditional diffusion equation) fractional 
Green function ( p. 21 ) reduces to the classic expression 



^'• t > = -^m^~m>- (2 - 22) 

Example 5. Wyss's fractional diffusion equation. 

The following example shows the effectiveness of the proposed method also 
for fractional integral equations. Let us consider the Wyss's type formulation 
of the diffusion equation || (for simplicity and comparison with the previous 
example — in one space dimension): 

u(x, t) = tp(x) + A 2 Dt a 9 U ^ ~ ■ (-°o<x<cx>, *>0); (2.23) 
lim u(x,t) = 0, u(x,0) =(p(x), (2.24) 

x— >±oo 

Applying the Fourier transform with respect to space variable x and the 
Laplace transform with respect to time t, we obtain: 

ff <A»> = £?& < 2 ' 25 > 

where U((3,p) is the Fourier-Laplace transform of u(x, t), is the Fourier trans- 
from parameter and p is the Laplace transform parameter. 

Inverting Laplace and Fourier transforms as it was done in Example ||, we 
obtain the solution of problem ( p. 23 ): 



/CO 
G(x-Z,tMOd£, (2.26) 
-oo 

1 f°° 

G(x,t) = - E al (-X 2 (3 2 t a ) cos Pxd/3. (2.27) 
vr Jo 



Let us evaluate integral ( |2.27j ). The Laplace transform of ( |2.27 ) and formula 
1.2(11) from |l{| produce 



oo 
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and the the inverse Laplace transform gives: 



G(x,t) = —J— / e p V f_1 exp(-|x|A-V /2 )*. (2.29) 

4A7T2 J 



Br 



Performing substitutions a = pt and z = \x\\~ l t~ p (p = a/2) and transfor- 
mation of the Bromwich contour (Br) to the Hankel contour (Ha), as it was 
done in a similar case by Mainardi 113], we obtain 



G(x,t) = / e*-™"-^- = -^rPM(z,p), z = ^- (2.30) 

K ' 2A 2vn J a l ~P 2A v Af V ' 

Ha 

where M(z,p) = W(—z, —p, 1 — p) is the Mainardi function ( |1.21|) . 

The last expression is identical to the expression, obtained by Mainardi |TI| 
in another way. 

We would like to note at this point, as in the previous example, that 
we have just evaluated the Fourier cosine-transform of the function U2((3) = 
E aA (-\ 2 p 2 t a ). 

For a = 1 the fractional Green's function ( [2.30 ) also reduces to the classic 



expression (gjp . The case of arbitrary number of space dimensions can be 



solved similarly. 

For a = 1 both generalizations (Nigmatullin's as well as Wyss's) of the 
diffusion problem give the standard diffusion problem, and the solutions reduce 
to the classic solution. However, it is obvious that the asymptotic behavior of 
( pl7D and ( ggg) for t -> 0, and t — > oo is different (see also discussion in on 



two different generalizations of the standard relaxation equation and discussion 



in [16] on two fractional models — one based on fractional derivatives and the 
other based on fractional integrals — for mechanical stress relaxation). 

This difference was caused by initial conditions of different types. The class 
of solutions is determined by the number and the type of initial conditions. 
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Chapter 3 

Sequential fractional 
differential equations 



3.1 The Laplace transform of a sequential 
fractional differential operator 



Let us consider the following initial value problem: 



Cty(t) = f(t); oDp-^t) 



t=o 



bk, k = 1, . . . , n 



n-l 



XMt) = a V^yit) + Y J Pk{t)aV' J t n - k y(t) + Pn (t)y(t), 
k=l 



(3.1) 
(3.2) 



3i 



1,2 



n ; 



0<ay<l, (i = 1,2,..., n). 



Fractional differential equation ( |3.ip is a sequential fractional differential 
equation, according to the terminology used by Miller and Ross To extend 
the Laplace transform method using the advantages of ( 1.15| ) for such equations 
with constant coefficients, we obtained the following formula: 

Laplace transform of the sequential fractional derivative QV^ m y(t). 

m—l 



poo "* 

/ e-P t V^f(t)dt=p^F(p) - p am -^' k -Dt m ' k f(t) , n , (3.3) 



a u t 



7-\^m — k 1 j~\^m— k — 1 



The particular case of (|3.3|) for f(t) m-times differentiable, a m = /j,, a*k = 1, 



(A; = 1, 2, . . . , m — 1) was obtained by Caputo IS, p. 41] much earlier. Taking 
«i = |U, 0^ = 1, (A; = 2, 3, . . . , m) leads under obvious assumptions to the classic 
formula ([2.1|). 
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3.2 Ordinary linear fractional 
differential equations 

In this section we give solutions of the "sequential" analogues of the "standard" 
linear ordinary fractional differential equations with constant coefficients. Of 
course, we must take appropriate initial conditions, also in terms of sequential 
fractional derivatives. 

Example 6. Let us consider the sequential analogue of Example |l|: 

oA° f aDlv(tj) + ay(t) = 



oD?- 1 ( D?y(t) 



t=o 



t=0 



&2, 



(3.4) 
(3.5) 



where < a < 1, < f3 < 1, a + f3 = 1/2. 

Our formula ( |3.3|) of the Laplace transform of the sequential fractional 
derivative allows us to utilize the initial conditions (|3.5|) . To use (|3~3|), we 
take a\ = a, «2 = P and m = 2. Therefore, o\ = a, o"2 = a + p. Then the 
Laplace transform (|3.3|) of equation (|3.4[) gives: 



y(p) 



(p«+0 + a)y(p)=p0&2 + &i, 

^ 1 



+ 6 



p a+ P + a p a+ P + a 



(3.6) 



(3.7) 



and after the Laplace inversion with the help of (1.15) we find the solution to 
the problem (Ot)-(O): 



y (t) = M Q -X+Aa(-a* a+/3 ) + ht a +P- l E a+/3t0l+ p(-at a +P). 



(3.8) 



For P = and a = 1/2 (and assuming, of course, 62 = 0), we can obtain from 
( p.8| ) the solution of Example |l|. 



Example 7. Let us now consider the following sequential analogue for the 

(3.9) 



equation, considered in Example [2]: 

oA a ( oD?v(tj) + oD q t y(t) = Ht) 



where < a < 1, < P < 1, a + P = Q > q. 
The Laplace transform (|3.3|) of equation 



gives: 

(p a+ ? +p q )Y{p) = H{p) +p?b 2 + b u 
h= oDr 1 ( D?y(t))\ t=o + Dr l y{t) 



(3.10) 



t=o 



t=0 



Writing Y(p) in the form 

p- q H(p) 



Y(p) 



+ b 



p 



p 



(3.11) 
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and finding the inverse Laplace transform with the help of ( |1.15| ) , we obtain the 
solution: 

y(t) = 6 2 t a - 1 ^ +/3 - 9 , a (-t a+ ^ 9 )+6i^-^ Q+/3 - g , Q+/3 (-t a+ ^' ? ) + 

+ f\t - T) a+ P- l E a+p _ q , a+p (-(i - T ) a+ P-«) h{r)dr. (3.12) 

It is easy to see that this solution contains solution of Example |asa special 
case. 

Example 8. Let us consider the following initial value problem for the sequen- 
tial fractional differential equation: 

D? 2 ( D^y(t)) - Xy(t) = h(t); (0 < ttl ,a a < 1) (3.13) 



oA" 2 oA ai y(*)_ = ki, oD^-'yit) t=Q = b 2 ; (3.14) 



1 wu, ■ in i i 

The Laplace transform (|3.3|) of equation ( 3.13| ) gives 



( p a 1+ a 2 _ X)Y(p) = p a2 b 2 + h, 

and after inversion using ( |1.15| ) we obtain the solution: 

y(t) = b 2 t^~ 1 E am {\t a )+b 1 t a - 1 E a , a {\t a ) + 

+ I \t-T) a - 1 E a<a (\(t-T) a )h(T)dT, (a = a 1 + a 2 ) (3.15) 
Jo 

Let us take a the same as in Example |. Using (|L1D , ( |L27| ) and flLlsD , it is 
easy to verify that ( 3.15[ ) is the solution of ( |3.13 ). It is also worthwhile to note 



that if 6i ± 0, b 2 ^ 0, then (gig) is not a solution of equation QDfy(t)—\y(t) 
h(t) from Example ^; also fl2.10| ) is not a solution of equation ( |3.13[ ). On 
the other hand, equations fl2.Sj ) and ( |3.13| ) are very close one to another: the 
fractional Green's function in both cases is G(t) = {t) a ~ l E a ^ a (\t a ) . 

3.3 Partial linear fractional 
differential equations 

Example 9. Let us consider Mainardi's |lf| initial value problem for the frac- 
tional diffusion-wave equation: 



D?u(x,t) = X 2 ^^- , (|x|<oo,t>0) (3.16) 

ox 



u(x,0) = f(x), (N<oo) (3.17) 

lim u(x, t) = 0, (*>0) (3.18) 

x±oo 



where < a < 1. 
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The Laplace transform of problem ( |3.16|)-(|3.18| ) using (3.3) for k = 2, a\ 
a — 1, a% = 1 (this gives Caputo's formula QlSfl ), i.e., 



C{ D?y(t)}=p a Y(p)-p a -ty(0) 



produces: 



p a u(x,p) — p a f(x) = A 



, d 2 u{x,p) 
dx 2 



(\x\ < oo) 



lim u(x,p) = 0, 



x±oo 



(t>0) 



(3.19) 

(3.20) 
(3.21) 



Applying now the Fourier exponential transform to equation (|3.2C| ) and 

i: 



utilizing boundary conditions (3.21), we obtain 



U(p,p) 



P 



p a + X 2 /3 2 

where U(/3,p) and F(j3) are the Fourier transforms of u(x,p) and f(x 



The inverse Laplace transform of the fraction 



(3.22) 



is E al (-\ 2 /3 2 t a ^ 



(E\.fj,( z ) is the Mittag-Leffler function in two parameters). Therefore, the inver- 
sion of the Fourier and the Laplace transform gives the solution in the following 
form: 



u(x,t) = J G(x-H,t)f(0dt, 



(3.23) 



G(x,t) 



7T 



1 



E a ^{-\ 2 f3 2 t a ) cos(px)df3 = ^-r p W(-z, -p, 1 - p), (3.24) 



where W(z, A, p) is the Wright function (|1.17 ). This solution is identical to the 
solution of the Wyss fractional (integral) diffusion equation ( |2.26| ). 
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Chapter 4 

Fractional Green's function 



4.1 Definition and some properties 

In this section we consider equation fl3.1| ) under homogeneous initial conditions 
b k = 0, (k = 1, . . . ,n), i.e. 

£ t y(t) = f(t); ^ fc ~Vt)[ =0 = 0, fc = l,...,n (4.1) 

n-l 
fe=l 

a^t — a-^t a-^t • • • a u t i a u t — a. u t o, u t ■ • • a^t > 

o-fc = 5I a i' ( fc = 1 > 2 ' • • • > n ); < «j < 1, (j = 1,2, . . . ,n). 

3=1 

The following definition is a " fractionalization" of the definition given in 

Definition. Function G(t, r) satisfying the following conditions 

a) T CtG(t, r) = for every r G (0, t); 

b) lim r ^i_ ( T Vl h ~ G{t, r)) = 5 k , n , k = 0, 1, . . . , n, 
(5k, n is Kronecker's delta); 

c) lim T ^ +Q ( T V^G(t,r)) = 0, k = 0, 1, . . . , n - 1 

T < « 

is called Green's function of equation (fOJ). 
Properties. 

Using (|1.28| ), it can be shown that y(t) = f G(t, r)/(r)dr is the solution 
of problem ( |4.1| ). 

Let us outline the proof of this statement. Evaluating oD^yit), oDl 2 y(t), 
. . . , oDf n y(t) using the rule ( 1.28Q and condition (b) of the definition of Green's 
function, we obtain: 



D?y{t) = oA Ql f G(t,r)f(r)dr 
Jo 
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= f T D?G{t,T)f{r)dT+ \my n rD^- l G{t,r)f{r) 

JO 

= f T D?G(t,T)f(r)dT (4.2) 
Jo 

D?y(t) = oA aa (oA ai )= oA" 2 /* rA Q1 G(t,r)/(r)dr 

JO 

= [* T D?( T D?G(t,T))f(T)dT 

Jo 

+ \im D^-\ T D^G(t, r))/(r) 

T^t — 

= [* T D?G(t,T)f(r)dT (4.3) 



oA*"" 1 ^) = oA Qn - 1 (oA CT "- 2 )= oA""" 1 f rD^- 2 G(t,r)f{r)dT 

Jo 

= f T D^( T D^G(t,T))f(r)dr 
Jo 

+ om r rl (,r 2 G(t,T))/(T) 

T^t — 

= /" T D?- l G{t,T)f{r)dT (4.4) 
Jo 

oA a "yW = oA a "(oA crn - 1 )= oA Qn f T D?-*G{t,T)f{T)dT 

Jo 

= f T D^( T D^G(t,r))f(r)dr 
Jo 

+ Km T D°»-\ T D?- l G{t,T))f{T) 

T^t — O 

= f T D^G(t,r)f(r)dr + f(t) (4.5) 
Jo 

Multiplying these equations by the corresponding coefficients and summa- 
rizing, we obtain 

oC t y(t) = f T C t G(t, T)f{j)dr + f(t) = f(t), (4.6) 
Jo 

which completes this proof. 

2. For fractional differential equations with constant coefficients G(t, r) = 
G(t—r). This is obvious because in such a case Green's function can be obtained 
by Laplace transform method. 
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The type (standard or sequential) of the equation is not important for deter- 
mining Green's function, because due to condition (b) in the Green's function 
definition all non- integral addends vanish. 

3. Appropriate derivatives of Green's function G(x, r) form a set of linearly 
independent solutions of a homogeneous (/(£) = 0) equation ( |3.1| ) (for a simple 
illustration, see Examples || and ^). 

Let us demonstrate this for the case of the linear fractional differential equa- 
tions with constant coefficients, which are the main subject for study in this 
work and for which we have G(t, r) = G(t — r). 

Let us take < A < a n . First, the function 

D?G(t) 



y\(t) 

is a solution of the corresponding homogeneous equation. Indeed, 



C t y x (t) = A( oA A G(i)) = oA A ( oAG(i)) = 0. 



A / 



(4.7) 



(4.8) 



We used here that o£* o-t^ = qD^ o£ t , which follows from condition (c) in the 
definition of the fractional Green's function. 
Second, 

-A-l 



D 



In fact, 

T->cr n — A— 1 



yx(t) 



i=0 







£}<X n — A— ] 



l yx(t) 



t=0 



1. 



t=o 



We first used here the relationship 



-A-l 



(cr„-A-l)+A 



G(t), 



(4.9) 

= 1. 
(4.10) 

(4.11) 



which follows from condition (c) of the definition of Green's function, and then 
condition (b). 

We see that having the fractional Green's function of equation (4.1), we 
can determine particular solutions of the corresponding homogeneous equation, 
which are necessary for satisfying inhomogeneous initial conditions. 

Therefore, solution of linear fractional differential equations with constant 
coefficients reduces to finding the fractional Green's function. After that, we 
can immediately write the solution of the inhomogeneous equation satisfying 
given inhomogeneous initial conditions. 

This solution has the form 



y(t) 



E b kMt) + G{t- T)f(r)dr, (4.12) 



i=0 



(4.13) 



Mt) = -D^~ ak G(t), V^-^ = a D^ aD?- 1 . . . a D? k+1 (4.14) 

Because of this, in the following sections we find some explicit expressions 
for fractional Green's functions, including a general linear fractional differential 
equation. 
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4.2 Fractional Green's function for the one-term 
fractional differential equation 

The fractional Green's function G\(t) for the one-term fractional-order differ- 
ential equation with constant coefficients 

a G D?y{t) = f(t), (4.15) 

where the derivative can be either "classic" (i.e., considered in the book by 
Oldham and Spanier) or "sequential" (Miller and Ross), is found by the inverse 
Laplace transform of the following expression: 

9i(p) = \ (4-16) 
ap a 

The inverse Laplace transform then gives 

Gi(t) = - S^- (4-17) 
a 1 (a) 



The solution of equation ( 4. 15| ) under homogeneous initial conditions is 



y (t) = — L- C - /(r) ^ r = - Dt a f(t). (4.18) 
yy ' aT(a) Jo (t-T) l ~ a a 1 J w v ' 

Using [28, lemma 3.3] , we can easily verify that expression Q4.18|) gives the 
solution of equation (|4. 15f) , if f[x) is continuous in [0, oo). 

4.3 Fractional Green's function for the two-term 
fractional differential equation 

The fractional Green's function G^i) for the two-term fractional-order differ- 
ential equation with constant coefficients 

a D?y(t)+by(t) = f(t), (4.19) 

where the derivative can be either "classic" (i.e., considered in the book by 
Oldham and Spanier) or "sequential" (Miller and Ross), is found by the inverse 
Laplace transform of the following expression: 

1 1 1 , 
92 (p) = = (4.20 

y yF> ap a + b a p a + ^ 

which leads to 

G 2 (t) = -t a - 1 E a ^{--n. (4.21) 
a a 

For instance, function (^(i) plays the key role in solution of Example [l] and S 
Taking in ([4~2l] ) b = and using the definition of the Mittag-Leffler function 



(|0|), we obtain Green's function G\(t) for the one-term equation. 
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4.4 Fractional Green's function for the three-term 
fractional differential equation 

The fractional Green's function Ga(t) for the three-term fractional-order differ- 
ential equation with constant coefficients 

a D$y(t) + b D?y{t) + c y(t) = f(t), (4.22) 

where the derivatives can be either "classic" (i.e., considered in the book by 
Oldham and Spanier) or "sequential" (Miller and Ross), is found by the inverse 
Laplace transform of the following expression: 

apP + bp a + c 
Assuming (3 > a, we can write (p) in the form 

1 cp- a 1 _ I VV-lV* f c \ k+1 p- ak ~ a 

93 ^ P ' ~ c avP- a + b CT>~ a ~c^' ' \a) (a h\ k+1 
cap +o l _w c k=Q \a/ [pp- a + i>) 

(4.24) 

The term-by-term inversion, based on the general expansion theorem for the 
Laplace transform given in []20|, §22], using ( [1.15 ) produces 

Gs(t) = It @ i ^ 1 H^(-^ 1 (4.25) 

k=0 

where E\„(z) is the Mittag-Leffler function in two parameters, 

We assume in this solution that a / 0, because otherwise we have the 
two-term equation ( 4.19Q . We can also assume c / 0, because for c = 

03 (p) = 1 = 

ap 13 + bp a apP~ a + b , 
and the Laplace inversion can be done in the same way as in the case of the 
two-term equation. 

Two special cases of equation ( [4.22 ) were considered by Bagley and Torvik 
fH (for = 2 and a = 3/2) and by Caputo @ (for (3 = 2 and < a < 1). It is 
easy to show that our solution ( [4.25 ) contains Caputo's solution as a particular 
case. 

Indeed, substituting (|4.26|) into Q4.25|) and changing the order of summation, 
we obtain: 

G3(t) ~ cg ( - 1} [a) g ( - 1} [a) Mj\ TW + k + 1) 



1 ^ yk(c\ k+1 (j + fc)l^+W-i-^ 

c^Ui^o ] \aj klj\T((3(j + k + l)-aj) { ' ) 

For [3 = 2 this expression is identical with the expression obtained by Caputo 
ftl formula (2.27)]. 
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4.5 Fractional Green's function for the four-term 
fractional differential equation 

The fractional Green's function G&(t) for the four-term fractional-order differ- 
ential equation with constant coefficients 

a D]y{t) + b Q D?y(t) + c D?y(t) + d y(t) = f(t), (4.28) 

where the derivatives can be, as in the previous section, either "classic" or "se- 
quential" , is found by the inverse Laplace transform of the following expression: 

04 (p) = nr ; ( 4 -29) 

K 1 ap~f + bp? + cp a + d v ; 

Assuming 7 > (3 > a, we can write g(p) in the form 
1 1 



34 ip) 



api + bpP l+ cp a +d 
ap" 1 + ' 



a^p-P 



p"< P + a 1 b g 1 cp a P+a 1 dp 13 
p^ 13 + a~ l b 

^ > (p7-^ + a-l 6 )m+l \ a P + a P 
\ ^ ( -i\m a P Y~W I ak~-/3m 

The term-by term inversion, based on the general expansion theorem for 
the Laplace transform given in |2(], §22], using ( 1.15| ) gives the final expression 
for the fractional Green's function for equation ( 4.28[ ): 

m=0 v 7 k=0 \ / 

(4.31) 

We assumed in this solution a / 0, because in the opposite case we have 
the three-term equation (4.22). We can also assume d 7^ 0, because in the case 
of d = 0, after writing 

54 (p) = 7 _ a ^.V a , ( 4 - 32 ) 

the Laplace inversion can be done in the same way as in the case of the three- 
term equation. 
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4.6 Fractional Green's function for the general 
linear fractional differential equation 

The above results can be essentially generalized. 

The fractional Green's function G n (t) for the n-term fractional-order differ- 
ential equation with constant coefficients 

a n D^y(t) + a n _! D^y{t) + . . . + a x D^y{t) + a D*>y(jt) = f(t), (4.33) 

where derivatives D a = qD" can be, as in the previous sections, either "clas- 
sic" or " sequential" , is found by the inverse Laplace transform of the following 
expression: 

9n(p) = a a " a 7T ( 4 -34) 

Let us assume (3 n > /3 n _i > ■ ■ ■ > fli > (3q and write g n (p) in the form: 
1 1 



E akpl 3 " 

CinP^ + dn-xp^- 



a n P 



1+ k=0 



g„-/3„-i | "n-1 

E / } n l lV n + i (continued) 



(m;feo,fci,... ,k n - 2 ) IJ ( ~ ) P 

ko + ki + ... + k n _2 = m i=Q 
k >0;... ,fc n _ 2 >0 



J/3i-i3 n -i)ki 



-y OO 



fl ™ m =0 k +k 1 + ...+k„_ 2 = 



^2 (-l) m (m;fco,fci,... ,fc n _ 2 ) (continued) 

) + fc 1 + ... + fc„_2=m 
fc >0;... .fc n _ 2 >0 

n-2 

-/?„_!+ E (ft-/3n-l)fe. 



where ((m; fco, fci, • • • , k n _2)) are the multinomial coefficients [22, chapter 24] 

The term-by term inversion, based on the general expansion theorem for 
the Laplace transform given in |2(], §22], using ( 1.15| ) gives the final expression 
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for the fractional Green's function for equation (4.33): 
1 00 (— l) m 

G n {t) = — ~ (m;k ,ki,... ,k n - 2 ) (continued) 

71 771=0 ' k + ki + ... + k n _2=m 

k >0;...,k n _ 2 >0 

JJ f 2) ^ ^-A-0"+/*.+E^o a (A-i-ft)*J-i (coutinued) 
i=o ^ a "^ 

ft,-ft,-i,m+EL(ft.-i-ft)*i V a n / v ; 
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